Photoabsorption in a plasma in a high magnetic field 
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Photo-absorption in fully ionized plasmas in high magnetic fields is re-examined, using the meth- 
ods of many-body quantum field theory. For frequencies in the immediate vicinity of the electron 
cyclotron resonance the rates we obtain disagree markedly from those in the literature. The new 
element in our work that causes most of the disagreement is the inclusion of the lowest order real 
part of the energy-shift of the resonant state, where, in effect, previous authors had included only 
the imaginary part. 
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Photoabsorption in a plasma in the presence of a strong 
magnetic field has been addressed many times in the lit- 
erature For the case of the dipole approximation 
to the absorption amplitude, and with the photon taken 
to be coupled only through the electron current, complete 
results can be found, for example, in ref. [lj. Recently 
there has been attention given to the effects of the proton 
current in a hydrogen plasma 

Gl-dH. This becomes im- 
portant in regions of extreme field and low temperature 
uj C p/T < 1, where u) cp is the proton cyclotron frequency; 
that is, when thermal photons have energy in or below 
the proton resonance region. We will follow others' ter- 
minology and call this the region of quantized protons. 
The present paper has three main objectives: 



a. To point out a complication that applies in either 
of the electron or proton resonance regions. Qualita- 
tively we describe it as follows: In the resonance region 
photoabsorption is best described in relation to resonant 
scattering. Then the usual free-free rate (times an easily 
calculated coefficient) becomes, in effect, a contribution 
to the imaginary part of the self energy of the propagator 
for the resonant state of the electron (or proton), to be 
added to the imaginary part coming from the collision- 
less resonance decay Of course, there is a real part to the 
self energy part as well, of the same order in the expan- 
sion parameter e 2 as the imaginary part. Since it is small 
we might have assumed that the real part is just an in- 
consequential shift of the resonance energy. However, it 
is energy dependent (and logarithmically singular at the 
position of the resonance). We find that including the 
real part greatly changes both the shape of the spectrum 
near the resonance and the integral over the resonance 
region. 



b. To exhibit an approach to the "quantized proton" 
problem that we believe is superior to those found in the 
literature, and to present a result for the "Gaunt factor" 
that is much simpler than that contained in the appendex 
of ref. 01. 



c. To address some collective effects that can be 
significant in domains of higher density. 

In discussing the issues that we shall raise, it is much 
more efficient, even for the recapture of single electron 
results or classical results, to begin from a quantum field 
theory formulation. This formulation is at the same time 
well adapted for addressing the emission, absorption and 
scattering processes that are involved, and for incorporat- 
ing the statistical mechanics of the plasma. We consider 
non-relativistic spinless electrons, of mass m annihilated 
by the field, V'e(r), and protons of mass M annihilated 
by the field ipp( r )- We choose Coulomb gauge for the 
electromagnetic field. Denoting the vector potential for 
the external field A(r), we define the current operator, 
as j(r, t) = j e ( r )t) + jp{ r ,t), where the currents for the 
respective fields {e,p} are, 
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V^HV + eA]xjj p + [zV + eA]V| V P . (1) 



We also define the number density operators n e , p (r) = 
^>|p(r)^ e ,p( r )- We divide the Hamiltonian into an un- 
perturbed part Ho, which includes all interactions with 
the externally applied magnetic field, 

Hq = — — j d 3 r\iV + el(r)}^ e • [-»V + eA(r)}yj e 
Zm J 



a Coulomb term, H c , which we divide into two pieces 
H c = H^+Hi B \ 

= 2 r 



HW = — 
8ir 



(d 3 r)(d 3 r')n e (r)-^—n e (r') 
|r — r'| 
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(d A r){d A r')n p {Y 



(d 3 r)(d 3 r')n e (r)- 



(3) 



(r'), (4) 
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and a radiation term, -ff ra d, which couples the external 
radiation to the matter. 



-Hrad 
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je "I" jp] ' -A-mcl 



(5) 



We choose the magnetic field B to be in the z di- 
rection. In this case it is convenient to define currents 
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± _ 



([je.p]x i i[j e .p]y) / 'a/2 which couple photon polar- 
ization vectors of the form (1, =pi, 0)/a/2. In the dipole 
approximation to photon absorption and emission we 
shall encounter the space integrals of the currents ([T]) 



J e , P = / d 3 r j e ,p(r) . 



(6) 



In a system governed by Hq each particle in the ensem- 
ble of electrons and protons moves independently and the 
second quantized formalism is unnecessary. For a single 
electron system, for example, the action of the operator 
J e on a state is the same as the action of the operator 
ell = e(p — eA) on the wave function for that state. Thus 
the familiar relations for an electron in a constant mag- 
netic field, [n+,II~] = eB = muj ce , [Ylf,H ] = ±uj ce Ilf , 
and [n+,n~] = -uj cp , [U^,H Q ] = Tw ce n±, translate 
into, 



[J^,Hq] — ^ui ce jf , [J^ , Hq] — ±uj cp J p 



(7) 



where Lo ce and uj cp are the respective electron and proton 
cyclotron frequencies. When we introduce the Heisenberg 
picture in the usual way, with H = Hq + H c + i? ra d; the 
equations ([7]) give, 

(£± Uce )j±(t) = [Jt(t),(H c (t) + H Tad (t))], 

T ^c P ) 4(t) = [J*(t), (H c (t) + H iad (t))} . (8) 

These are the key equations for our application. Turning 
to the commutators of the J's with the Coulomb Hamil- 
tonian, it is more transparent if we use the 3D vector 
representation and express the J's in terms of infinites- 
imal translation operators for electron coordinates T e , 
and for proton coordinates, T e , 

Je = — T e + e 2 [ d 3 rn e (r)A(r) , 
m J 

J p = |^T p + e 2 J d 3 r n p (r) A(r), (9) 

where, 

[T e , p ,^e,p(r)] = VV> e ,p(r). (10) 
Noting that [n e p (r)A(r), H c ] — we find first that, 

[J e , p ,HW]=0, (11) 



since the e-e and p-p interaction terms are separately in- 
variant under the separate translations of either the elec- 
tron or proton coordinate. However the e-p interaction 
He is invariant only under simultaneous translations 
for the electrons and protons, 



[T e + T p , H c 



M 



0. 



(12) 
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The separate commutators needed in ©are thus, 

M 

[J e ,Hc\ = [J p ,Hc] 

m 

J(d 3 r)(d 3 r')V r n e (r,t)y^n p (r',t). (13) 
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FORMAL EXPRESSIONS FOR THE 
PHOTO-ABSORPTION RATE. 

Since the issues raised in the present paper pertain 
largely to the modes polarized perpendicularly to the 
magnetic field, we address only the absorption of these 
modes. We begin from the expression that the field the- 
oretic formulation of statistical mechanics provides for 
the photo-absorption rate 7a(w)± in the medium, where 
± refers to the photon modes with polarization vector 
(l,±t,0)>/2 



7a(w) ± = 2^kql d3xd3 y dt 
x(j±(x,i)j T (y,0)) 



e zw* e -iq-(x-y) 



(14) 



The Heisenberg operators ^ (x, t) are the relevant 
combinations of the total electromagnetic current oper- 
ator, j ± = ([j e + j p ]i ± i[j e + j p ] 2 )/V2. The brackets, 
() stand for thermal average, (A) = Tr[e-' 3H A]/Tr[e-' 3H ] 
where f3 — T _1 . Taking the dipole limit, q — > 0, and 
using the definition = J d 3 rj e t p (r) we obtain, after 
an integration by parts, 

7^H ± = £;Jdte^^(§±u; ce )(Jt(t)j*(0,0)) 
+ijJdt e<" ^ (I T c cp ) (J* (t) r (0, 0)) ■ (15) 

Now we substitute ([5]), letting A lad — since (at the mo- 
ment) we are considering a process in which no photons 
are involved except the one being absorbed. 

TaM* = & I dt e** ( [{ ^ + ^ } , H c (t)] (0)) 

= _i r d t e ™ ( \t'im + sm. \ , h c (o)} w> • 

(16) 

The second form follows from translational invariancc 
(to shift the space integral from the first to the sec- 
ond current), space inversion (introducing a (-) sign) and 
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time translational invariance (a displacement —t) in the 
thermal average factor, followed by complex conjugation. 
Performing the same steps once again we obtain, 



7A(«) ± = ^/Ae^<[(^ + ^£) tJ ff e (t) 

*[(tS + S^Wo)]>. 

Using © and (jlll) we can rewrite (fTTHm a simpler form 1 

1 n2 



(17) 



7aM 



[Vol.1 VMfwi 



2w[Vol.] \M(u ± ui C p) m(uj =p w ce ) 



eft e l 



T±(t),H c (t) Tf(0),ff c (0) ) 



(18) 



or explicitly, using (Ti"2"]) 
7a(w) ± 



m(o;±a; C p)(a;=pa; ce ) 



x /(d 3 n)(d 3 r 2 )(rf 3 r 3 ) 



J-r fdte* 

r:5 J 



|ri-r 2 | |r 3 

x (d±n e (r 1 ,t)n p (r 2 ,t)d T n e (r 3 ,0)n p (0,0)) . 

(19) 

where we have defined, d± = (d/dx± id/dy)/y/2, and 
simplified the prefactor using Mu> c 



moj c 



and 



(to + M)/M sa 1. This prefactor in (fT^j) agrees with 
that found in ref. (5j, except for the damping terms that 
matter only very near the resonance in the results of these 
authors. This subject will be discussed at length in sees. 
4 and 5. 

To recapitulate, (fT9"f gives the exact rates in the dipole 
limit. All of the effects of Coulomb interactions in the 
medium are included. Of course, the four-density cor- 
relator must calculated in an approximation. There are 
different domains in which it makes sense to take dif- 
ferent approaches. For example, at high temperatures 
and low densities we can use the result in the absence of 
Coulomb coupling (beyond that already exhibited explic- 
itly in (1191) ). For low temperatures and low densities we 
need to include full Coulomb wave-functions in the elec- 
tron's interactions with a single ion. For high densities 
collective effects become interesting, and "ring approxi- 
mation" sums required (at the least). 

BORN APPROXIMATION IN TWO REGIONS 

The four-point density correlator depends on both 
Coulomb and magnetic interactions. By "Born approxi- 



1 The volume, Vol. , enters or leaves our formulae depending on 
whether or not we use, e.g., the current evaluated at an arbitrary 
point in space, j'J^O, 0), as in 1171 or the space integral of the 
current which, as it occurs in H18I I. is just the current at the point 
times the volume, in view of the translational invariance of the 
thermal average factor in dlTI . 



mation" we mean that the Coulomb couplings are turned 
off in calculating the bracket in (fl"9]l . In this case the cor- 
relator factors into an electron part and a proton part, 
separately translationally invariant, 

(d±n e (r 1} t)n p (r 2 ,t)d T n e (r 3 ,0)n p (0, 0)) 
= (S±n o (ri,t)a F n e (r 3 ,0))<Ti 1) (ri,t)Tj 1 ,(0,0)> . (20) 

Using (|4"Tj) and taking Fourier transforms in (fT9"]l gives 

taH^-Gt) 3 

W7r\47r/ 



m(u ± w cp )(w =F ^ce) 



F(lo) 



(21) 



where, 

F(u>) = J dwi d 3 k 
and 

A e , p (k, u) = J 



kl 



(fc 2 ) 2 



A e (k,w-wi) A p (k,wi) 



(22) 



d 4 I e* x e « w '(ne, P (x,t)n e , P (O,0)) . (23) 



to be evaluated in the absence of Coulomb interactions. 
We work out (|22p in two domains; ui cp « T < 1.5w ce 
and T « 1.5w cp , the first being a domain in which the 
electron occupancy is mainly confined to the ground and 
first excited Landau levels, with the magnetic effects on 
the proton being negligible. The second domain is one 
in which the electrons are strongly confined to the lowest 
level and the protons are mainly confined to the ground 
and first excited Landau levels. 



Region: lj cp << T < l.5ui ce 

In this domain the free proton correlator is simply, 

A p (k,w) = 27rr4 0) c5(w) . (24) 

There are collective effects that modify this correlator to 
which we return later. We introduce the variables, 



Si ~ V 2 



t — HIE j£ ZL 

S° V 2 \ fen 2m 

£-1 



m/3 f u; 
2 



2m 



c 



2mto c 



2m j ' 

(25) 



Following the rules given in the appendix we have calcu- 
lated the terms in the electron correlator that come from 
the first two Landau levels only, 

A e (k\\,k x ,u>) = (1 - e-^^^/^ivm^nf^k^e-i 
Cexp(-C?) + [1 + (1 - C) 2 e-^"]exp(-£ 2 ) 



+e-^«Cexp(-£i 1 ) 



(26) 
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Substituting (|26|) and (p4|) into (|22|) . introducing a vari- 
able s = fc^/fcy, and then doing the k\\ integration, we 
obtain, 



F( W ) = H A /2^8 7r 3/2 [n (o) ] 2 A) 



(27) 



where A has exactly the same meaning as in ref . [B[ , and 
is given by, 

A=f(l-e-^)e^/ 2 / o °d SlT ^ F 
x K (puy/.25 + S /u C eP) (l + £3s±d*-P»c 



(s+1) 2 



+ (1 + S )-l e -^ee/2 [K (P\U - W ce |V-25 + S / Wce /?) 

+#0 + w ce | y.25 + s/wce/3)] ] . (28) 

This answer is identical to the results of Pavlov and 
Panov [H, as corrected by Potekhin and Chabrier in 
eqn. (44) of ref. [j| , when the latter is expanded in powers 
of exp[— jScJqb], and only the zeroth and first order terms 
retained. 



Note the very close resemblance to (|28|) . even though 
([30)1 is to be used in a domain of temperature 1000 times 
smaller, at a given magnetic field. Of course when we 
took the two lowest Landau states for the case of quan- 
tized protons, rather than for electrons, the energy dif- 
ference u! ce in (|28p is replaced by cj cp in (|3"0"j) . But note 
that it is still the electronic parameter uj ce that enters 
the y/.25 + 2s/u ce factor in the arguments of the Bcsscl 
functions, , but with a coefficient that is different by a 
factor of 2. When the temperature is so low that we have 
uj ce f3 — 2000, then the integral of the first Kq function in 
(l30l) becomes rather large, since convergence for large s 
comes from the cutoff supplied by the K function. 

Potekhin and Chabrier [5] have given formulae which 
should exactly agree with (f30|) ■ However, they went quite 
a different route to obtain these formulae and end up with 
complex expressions that we have not been able to cast 
into our form. However, we obtain a plot similar that of 
the lower panel in their fig. 6. showing the peak in the 
scattering rate for the anti-resonant polarization, when 
plotted over the proton resonance region. 



Region: T < 1.5uj c 



We choose the region to extend to 1.5u> cp in order 
to capture the resonance behavior, while keeping the 
electrons strongly confined to the lowest Landau level. 
Therefore we take only the term with unity in the final 
factor in ([2"B)) for the electronic correlator.. The contribu- 
tion of the first two Landau levels to the proton correlator 
is given by taking m — > M, and lu cc — > uj cp in (|26[) as it 
stands. Doing the u>i integral in (|2"TT) . discarding terms 
of relative order m/M , and setting exp(— f3uj ce ) = 0, ap- 
propriate to the temperature regime, we obtain, 



F(u) = (1 



xe 



-2C 



x exp 



exp 



^)v^fl7r[ni 0) ] 2 / 



d*k 

TfcTT 



51 



— /3w c 



2k' 2 



exp 



(fcjj+fci) 2 
(i-0 2 ) 

(a;— uj cv ) 2 mj3 
2fc? 



+Ce- 



exp 



(29) 



In the calculation we replaced the reduced mass by m 
in several places. Then the only place the proton mass 
enters is through uj cp . Note that ( = k\/ (miu ce ) = 
kj_/(Mu) cp ). Doing the kn integration we obtain (|27|) 
with A replaced by A', where 



A' 



(1 



+(1 

-(n 



b(j8|t 



2{s+\) 2E 

2s/u; ce pJ 



,1/25 



■s)- x e-^' 2 KQ[P\w 



1^.25 + 25/^/3 



(30) 



ELECTRON RESONANCE REGION. 

We now address aspects of behavior in the resonant 
region that we believe are not adequately treated in the 
present literature. Specializing to the case of our first 
domain uj cp « T < 1.5w ce appropriate to the electron 
resonance region, ([2Tj) becomes 



7a(w) 



64tt 4 o 



m(cj - LJ ce ) 



(31) 



an expression that doesn't exist at the resonance fre- 
quency; we must turn to the u) — > lu — it prescription for 
the definition. The way that this works is that we first 
recognize (I31[) for the non-resonant case as the imaginary 
part of an amplitude defined by the graph of fig 2. 




FIG. 1: The basic 7 + e — > 7 + e graph, the imaginary part 
of which gives the resonant part of the phot oabsorpt ion rate. 
The heavy lines indicate the first excited Landau level. 
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Accordingly, we replace (l3lj) by 



7a M 



64tt 4 o; 



Im 



m(oj 



G{u)-ie) . (32) 



where -F(w) = IniG(a; — ie). When we move into the 
resonance region the denominator factor (uj — uj ce ) must 
have an imaginary part as well. In our original defini- 
tion of the problem in which the only photons are the 
external ones that coupled to the currents in (TT41 this 
imaginary part originates in a self energy insertion, 
in the inverse propagator of the resonance. The imagi- 
nary part, ImE (//) comes from our earlier calculation, 
but redescribed (with the appropriate multiplying coef- 
ficients) as the collisional deexcitation rate, Vff for the 
resonant state. Of course, the resonance inverse propa- 
gator also has an imaginary part coming from the inter- 
mediate state in which a photon reappears, defining the 
radiative (or natural) width v re — (2/3)e 2 uj ce uim~ 1 in 
the absence of Coulomb collisions. We denote the sum of 
the two imaginary parts, Imp^) + E^], by v e , where 
Ve = Vff + v re . 

It is conceptually incorrect simply to replace (ui — 
uj ce )~ 2 in (pJ2"j) by [(ui — w ce ) 2 + v 2 ]^ 1 . It also can lead to 
very incorrect numerical results, depending on the mag- 
nitude of Re[G]. Instead (p?2"j) should be replaced by 



7aM 



2 (0) 



Tm[u - uj ce - E« - E^]" 1 , (33) 



where v e = ImpM + S^']. The prefactor in (J33]) is 
determined by expanding to first order in ImE^) = Vff 
and comparing to (|31j) . with the identification, 



4 ImS(W>(w) _ a 2 F(uj) 



QAir A mu) ce n 



(0) 



47r 2 ma; ce ni ^ 



(34) 



A graphical representation of (|33|) is shown in fig. 2. 

The result ®, where we use E^ + £W) = i{v ff + 
Pre), gives exactly the usual results [5] for the resonant 
region. Now, however, we can ask about the effects of 
YfJfi and £M If these real parts of are small and not 
too energy dependent, then they provide small shifts in 
the resonance energy that make no difference either to 
the total rate or to the shape of the spectrum. But it 
turns out that E^)(w) has a logarithmic singularity in 
its real part at oj = uj ce that is of the same nature as 
that in its imaginary part, and we will see that it can no 
longer be dismissed. 

In calculating ReE^), we can avoid introducing a for- 
malism with propagators for the resonant state by di- 
rectly writing 



S (//) 



a 



irmu> c 



d A k 



[n e (k,u;)] 



(fc 2 ) 2 J L(l-e-^ w ) 



(35) 




FIG. 2: A typical graph that enters the photoabsorption 
calculation in the resonance region. The incoming photon is 
absorbed on an electron in the lowest Landau level (exciting 
the electron to the next Landau level). The initial state is 
restored at the right hand side of the diagram. In between we 
have the propagator for the resonance, as given in (|34[) . with 
self energy insertions corresponding both to the "free-free" 
scattering from protons and from free decay into photons as 
well. The photoabsorption rate is found from the imaginary 
part of the sum of all such graphs. 



where we have used (]22j) and ([24]) . and have extended 
A(k, a;) of (|2"2"j) to the complex plane with, 



A(k,w) 



Im[n e (k,h;)] 

(l- e -/J«) ' 



(36) 



The imaginary part of (|35|) reproduces the above results 
for Vff. The function II e (k, ui) is the Fourier transform 
of the retarded commutator (see ref.[l2j], sec. 33), 



II e (r,i) = ([n e (r,t),rae(O,O)]>0(t). 



(37) 



which supplies the extension to the complex plane that 
has the correct analytic properties. 

To construct the real part of II we extend the variable 
set (ESJ), 



f ± _ ImP ( uj_ , J^_\ a± _ / m/3 
^0 — V 2 \k n ^ 2m J ' SI — Y 2 



x 2m/ i 



C± _ m/3 
S2 — A/ 2 



I 2m J ' ^ 2muj c 



and we define the usual plasma function, 



$(0 = 2e- f3 / dy< 
Jo 



(38) 



(39) 



Then we find, 



Repl e (k,w)] = 



/ m/3 n e e ^ (1 — e 
V2ku 



x[($(e Q -)-^ + ))(i+ e -^(i-c) 2 ) 



(40) 
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With the real part of of (|35|) determined from 

(|40]) and the imaginary part from (|36|) the absorption 
rate is calculated from from ([55]) . In fig. 3 we show the 
difference that inclusion of the real part of the resonance 
self energy term can make in a region near the resonance 
peak., for the case of parameters, p = 10<?c~ 3 , u> ec — 1 
KeV,T = .2KcV. Plotted are the opacities in a region 
l.OlKcV < uj < l.lKeV, coming from ([33]) using the real 
part of T,^ that comes from (|40|) and for the imaginary 
part simply Vff as computed in section 3. Neither the 
real nor the imaginary parts of £( re ) make a significant 
contribution in this region. 
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In the first term the photon is absorbed on the electron; 
since we are totally out of the resonance region we do 
not need to build it into the propagator for the electron 
resonance, as above, although there would be no harm in 
doing so. We evaluate it as is, denoting the contribution 
as 7(1), 



7(i) 



OJ,, 



7r(47r) 3 wM 2 cp 



(42) 



The cross term, that is, the interference between the am- 
plitude in which a photon is absorbed on the electron 
and that in which the photon is absorbed on the proton, 
is trickier. It is not a constant times the imaginary part 
of a self-energy in a resonance propagator; we evaluate it 
almost as it stands. We recognize that it is a piece of an 
imaginary part of a function of something with four legs, 
but save ourselves the chore of further formal definitions, 
by noting that the function F(ui) defined in (|2"21 is the 
imaginary part of an analytic function, G(uj), the real 
part of which we have calculated, but which is insignif- 
icant in what follows. The contribution of the second 
(cross) term in (|4Tj) becomes, 



FIG. 3: Comparison of the absorption rates, ja, very near 
the electron resonance with and without inclusion of the real 
part of the self energy part in the resonance propagator. The 
dashed curve is the conventional model with no real part. The 
resonance energy is .2 KeV and the plot is from 1% below the 
resonance energy to one percent above. 



PROTON RESONANCE REGION 

With these results of the previous section in mind, we 
return to the case in which both proton and electron cur- 
rents contribute, where the interference term is impor- 
tant. We will consider the ^ polarization state, which 
resonates with the proton at energy uj cp , but never with 
the electron, and stay in the energy region lo << u ce 
throughout. We go back to (|2"T|) . rewritten to separate 
the proton and electron current as in (TIB)) . 



, * _ e 6 F{uj) / 1 + 1 \ 2 

w(47r) 3 V M (u) — UJrri) mild + OJ r p) ' 



e e F(oj) 
jM 2 (4tt) 



(47r) 3 \M(u- Lo cp ) m(u + uj ce ) 
uj^p + 2w cp 1 (u; - lo cp )~ x + (w - u cp ) 



(41) 



where we have used muj ec — Muj PCi and u << Lo ce to get 
the second form. 



7(2) M 



2e b 

(47r) 3 ww cp M : 
2e 6 



rlmf 



G(u>) 



(w - u cp )F(uj) 



(4ir) 3 uiuj cp M 2 (to — oj cp ) 2 + u p (ui) 2 ' 



(43) 



where v is just the imaginary part of the proton- 
resonance inverse propagator under the influence both 
of the magnetic field and the Coulomb scattering, v p = 
v rp + Vff with (2/3)e 2 M" 1 uJOJcp. In the second 

line of (|4"3"]) we have set the real part of G equal to zero 
after all, after calculating it approximately and finding it 
inconsequential. (It was important to have it included in 
principle, since our calculation should always be that of 
calculating the imaginary part of an analytic function.) 

As we found in the last section the free-free width term 
in the denominator Vff can be expressed in terms of the 
function F of (l22las 



a 2 F(uj) 



ATt 2 mio ce ne 
where for this case F(ui) is evaluated as 

FH = ^V2^ 3/2 k 0) ] 2 A', 



(44) 



(45) 



where A' is evaluated in (|30|) . 

The third term, in pT]) leads to a contribution ja (w) (3) 
where the factor [u> — oj cp )~ 2 needs to be fitted into the 
propagator for the proton resonant state, exactly in the 
fashion used in the pure electron problem to obtain (|33[) 



7 



7(3) 



2 (0) 

e n e u), 



cp 



ujM 



lva[u) — u> c 



(46) 



We can now add the three contributions, replacing (|4"Tj) 

by, 1A -> 7(i) +7(2) +7(3), 



ujuj cp M 



(UJ -UJ cp ) 2 +1/2) 



(47) 



Again we appear to disagree with the results of ref . [5| . 
According to eq. 53 of that paper, the free-free part of 
the damping term corresponding to v p as it occurs in the 
denominator of (|4"T)) . above, carries an additional factor 
of m/M compared to our expression. Also the second 
term in the numerator of (|47[) is missing in ref. 



COLLECTIVE EFFECTS 

For problems involving Coulomb forces in a plasma, 
screening is the leading collective effect. Indeed the au- 
thors of ref. [7[, and previous authors as well, in effect 
replace the factor [fc 2 ]~ 2 in (|2"2"|) by the Fourier transform 
of a screened potential, [fc 2 + k%]~ 2 , where kd is the 
usual screening parameter. This replacement is already 
known to be incorrect in the B = case; the correct 
static screening correction for photo-absorption in a hy- 
drogen plasma has been shown 14{ instead to be the 
replacement, 



-> fc" 



,1 



<D 



/(2fc 2 ) 



(48) 



We expect the same result for the magnetic case; in 
accord with Sitenko's conclusions 15], but we shall 



nonetheless look at the matter in some detail, as there 
are non-static corrections that may be significant in some 
regions. We return to (|22|) but re-express the right hand 
side using (f3"Tf to obtain, 



F(u>) = J dux d 3 k 



Im 11^ (k.w — oji) Im lip (k.ui ) 



(i- 



)) (x-e-f^i) 



(49) 



where the II' are the functions defined by (13"T|) . but now 
calculated, in approximation, in the presence of Coulomb 
forces, whereas the the functions, II are the functions 
defined in (|37|)in the absence of Coulomb forces. In the 
ring approximation, we have, 



n' e (k,a,) = 

n;(k,^) = ¥ 



[fc 2 +47re 2 n p (k,w)]II e (k,u;) 



fc 2 +47re 2 n e (k,a;) +47re 2 n p (k,tj) ' 
[fc 2 + 47re 2 n e (k, uj)}U p (k, u) 
-47re 2 n p (k,tj) +47re 2 n e (k,a;) ' 



(50) 



We can we describe this construction as the expression of 
the complete polarization parts II' e ,II p from the proper 
polarization parts, II e ,IIp, then setting the proper parts 
to their values in the absence of Coulomb interactions. 
The best systematic derivation that we know of for these 
relations in a multicomponent classical plasma is in ref. 
13], cq. 2.110. 

The relations hold in the presence of quantum effects 
as well. We have obtained tractable expressions from 
substituting the results of (|50|) into (|49|) only for the first 
case of section 3, uj cp « T < 1.5oj ce , where we can 
neglect the magnetic interactions of the protons. For this 
case the large proton mass leads to the imaginary part of 
the proton correlator Im[D p ] in ([21]) being concentrated 
at very small values of uix << T, for relevant values of 
fc « y/mT. Thus in (|21[) we can set u> — lo\ = lj, and 
(1 — exp[— flujx])^ 1 — (/?wi) _1 , and use the dispersion 
relation to do the cox integral, 



ir- 1 ^- 1 J dw 1 lmll' (k^O/wi =^- 1 ReI^,(k,0) 



_ 47r e 2 »(°)(l + B ^/fc 2 ) 

(l + K2/fc2 +K 2/ fe 2) ! 

where we have simplified by setting 

47re 2 n e! p(k, 0) w 47re 2 n e!P (0, 0) 



(51) 



(52) 



Here the k 2 p are the contributions of the individual 
species to the squared Debye wave number k 2 d = k 2 + n p 

and n p — k 2 — 47r/3e 2 ni°' ) . We have also rewritten the 
multiplying factor of K p in (|43[) terms of the average pro- 
ton density, n p ^ — and the temperature, ■ For 
the static response function (lj = 0) this is a good ap- 
proximation for all cases under consideration. 

Putting these steps into (|2"Tj) , using the ring approxi- 
mation (I50p for the electron polarization, and noting that 
Hp(k,uj w (i^ 1 ) w gives 



F(oj) = - 



Jd 3 k 



Im n< 0) (k,w) 



|i+nr' (k,^)/fc 2 | 2 



(X+nj/k 2 ) 



(X + Kl/W+Ki/W) ■ 



(53) 



We can compare f|53[) for the case B — with re- 
sults in the literature giving the effects of Coulomb cor- 
relations on the photo-absorption rate. Taking the ap- 
propriate limits of the prefactor in (|19[) , substituting 
the B = form for Imlli ' ) /fc 2 , and defining e(fc,w) 
1 + Ili ^, w)/fc 2 , and including the longitudinally po- 
larized modes gives back exactly eq. (3) in the paper 
by Iglesias and Rose [3], and essentially the results of 
ref. [17| as well. 2 The corrections from the |e| 2 in the 



2 There is one discrepancy; in both of these references the dielec- 
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denominator and from the ionic correlator are actually 
relatively small in domains of density and temperature 
in which the plasma is weakly coupled, that is to say, in 
regions in which we can calculate at all. We note that 
the last factor on the right hand side of (|53|) provides the 
screening factor that we quoted at the beginning of this 
section. 



DISCUSSION 

We summarize our differences from previous authors 
in three different parameter regions: 

1. Region of non-quantized protons. Here we re- 
capture the results of other authors for the "logarithm" , 
or "Gaunt factor", A. But there is another effect, at 
the same level of approximation, that can change the 
shape of the resonant peak by a lot, namely the rapid 
energy dependence of the real part of the self energy in 
the propagator for the resonance. Because of this energy 
dependence the real part is not merely a small adjust- 
ment to the resonant energy parameter, and can make 
large modifications to both the spectrum shape and the 
integral over the resonance region. We believe that the 
effect will be important in the proton-resonance region 
as well, but we have not calculated examples. 

2. Region very close to the proton resonance. 

While we agree with the authors of ref. [5;] in the way 
that the free resonance decay partial width parameter 
v rp enters the total width parameter v pi we appear to 
disagree with these authors on the contribution of the 
free-free parameter Vff itself to this width. 

3. Collective effects We have included the basic 
ionic screening effects which are somewhat different from 
those of other authors, and which, in contrast, reduce to 
well known collective corrections to photoabsorption in 
the limit of no magnetic field. 

APPENDIX 

The electron field operator ip(r, t) is built in cylindrical 
coordinates from the states of the theory, in the usual 



trie function e, where it occurs in the denominator, is taken to be 
that of a classical plasma. The numerators, which are not in these 
works identified as the imaginary part of the dielectric functions, 
require the quantum treatment in order to avoid an ultraviolet 
divergence (or the introduction on an arbitrary logarithm). In 
any application the full quantum form should probably be used 
in the denominator as well. 



way, 

where p is the momentum in the direction of the field 
&n,s{p) is the annihilation operator for the indicated 
mode, and e p = p 2 /2m. To calculate the correlator 
needed in (|23j) . 

A e (x,i) = (^,^,^(0,0)^(0,0)) , (55) 

we first write 

^(r,^ f (O,0) = -^(O.O^M) +C(r,t) (56) 

where the anticommutator function C is a c-number. In 
the present work we are considering only non-degenerate 
electrons and therefore we discard the first term on the 
RHS of (1561) . since its contribution is of higher order in 
the fugacity of the electrons, e" Me , i.e., a correction for 
Fermi statistics. (Note that the function C itself would 
be unchanged for the case of Bose statistics, while the 
discarded term would be of the other sign.) 

Before proceeding further we note that since the func- 
tions it n ,s(i"j_) for n ^ s vanish at rj_), and one elec- 
tron field is evaluated at rj_) in each expression that we 
encounter, only terms with n = s will enter, and the 
azimuthal angle cf> will not appear in any expressions. 
Henceforth we use radial functions labeled with the in- 
dex n alone; u n (p) = u n ^ n (p). 

Explicitly, the function C is now given by 

C(\r x \,z,t) = J2 I ^e^e-^^+^VCKIKCO) 
n *' 77 

(57) 

For the remainder of the evaluation of (|55l) and its Fourier 
transform we need the thermal expection value, 

([a n (p)^[a n (p')]) ^ 6 n ^S p , p :e^^e-^/ 2m+n ^ , (58) 

where to determine the fugacity exp[/3/i e ] to be used in 
([55]) we calculate the electron density (for convenience at 
r = 0) as 

n c =e^Y\u n {0)\ 2 e- n ^ ± e -^/(2 m ) (59) 
Using |m„(0)| 2 = TOW c /(27r) we obtain, 

_3/2o3/2 (0) al/2 

= (1 - e-""") - 2 " e g (60) 



Putting together (|56j) . (l57l) . and (|6T1) we obtain 
A e (k,w) = e^" J d 2 r ± dze 



oo 

X 

3=0 



E / ^e^ % (|rx|)%(0)C(rx^,t) (61) 



For examination of behavior in the region of the main 
cyclotron resonance it suffices to include only the terms 
j = 0, 1 in the sum in (|61|) and terms with n = 0, 1 in the 
sum in (|57[) . Inserting 



Ul (p) = ^= ^e"~ ' 4 (62) 

where 7 = oj c m/2, and performing the Fourier transforms 
yields the result ([2"rJ]l. 
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